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Theory  of  electron-avalanche  breakdown  in  solids 


★  ★ 

M.  Sparks,  D.  L.  Mills,  R.  Warren,  I.  Holstein, 

★  ★  ★ 

A.  A.  Maradudin,  L.  J.  Sham,  E.  Loh,  Jr.,  D.  F.  King 

Pacific-Sierra  Research  Corporation 

Santa  Monica,  California  90404 


Electron-avalanche  breakdown  in  solids  is  explained  by  a  theory 
that  is  predictive  and  agrees  with  experimental  results  for  the  magni¬ 
tude  of  the  breakdown  field  and  its  temperature  dependence,  pulse- 
duration  dependence,  material-to-material  variation,  and  wavelength 
dependence  for  X  good  agreement  between  experiment  and 


Jtyie  good  agreement  between  experiment  and 


theory  with  no  parameters  adjusted  is  obtained  by  using  improved  mag¬ 
nitudes  and  energy  dependences  of  the  electron-phonon  relaxation  fre¬ 
quencies.  The  contributions  of  both  optical  and  acoustical  phonons 
to  electron  loss  and  energy-space  diffusion  must  be  included.  The 
breakdown  field  tjg^is  calculated  by  solving  an  eigenvalue  equation 
obtained  from  the  diffusion  transport  equation.  Simple  models  and 
interpretations  of  the  diffusion  equation  afford  physical  insight 
into  breakdown  and  render  the  breakdown  conditions  predictable.  Pre¬ 
liminary  results  indicate  that  the  diffusion  approximation  fails  for 
wavelengths  considerably  shorter  than  1  um,  where  multi  photon  absorp¬ 


tion  must  also  be  considered. 


y V !  1  C*\  P  O 


I.  INTRODUCTION 


Damage  of  dielectric  materials  by  high  electric  fields,  which  has 

been  the  subject  of  numerous  experimental  and  theoretical  investiga- 
1  2 

tions,  ’  has  received  renewed  attention  in  recent  years  as  a  result 

2b 

of  laser-damage  experiments.  Electrical  breakdown  of  dielectrics 
by  lightning  was  one  of  the  first  known  electrical  phenomena;  labora¬ 
tory  experiments  on  electrical  breakdown  in  glass  were  initially  con¬ 
ducted  in  1799. 1  Despite  many  previous  attempts^-17  to  develop  a 
viable  model  of  laser-induced  damage  to  alkali  halides,  there  has  yet 
to  appear  a  theory  capable  of  simultaneously  explaining  the  magnitude 
of  the  breakdown  field  as  a  function  of  temperature,  wavelength,  pulse 
duration,  and  variation  from  material  to  material. 

Predictions  of  the  electron-avalanche  theory  and  a  comparison 
of  theoretical  results  with  data  available  at  the  time  were  given  by 

Sparks J 7  A  recent  account  of  theoretical  efforts  by  Soviet  analysts 

12 

was  given  by  Gorshkov  and  coworkers,  and  new  data  taken  at  the 

Lebedev  Institute,  along  with  further  discussion  of  the  theory,  were 

1 3 

presented  by  Manenkov.  Difficulties  with  the  electron-avalanche 

14 

theory  led  Schmid  and  coworkers  to  explore  such  alternative  mech¬ 
anisms  as  multiphoton  absorption  combined  with  a  possible  rapid  trans¬ 
fer  of  electron  energy  directly  to  the  lattice  by  phonon  emission— 

15 

following  a  suggestion  of  Hellwarth  that  direct  and  rapid  transfer 

of  energy  to  the  lattice,  followed  by  its  fracture,  may  account  for 

the  data.  Preliminary  results  of  the  present  investigation  were  sum- 

1 8 

marized  briefly  elsewhere. 
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The  general  features  of  the  present  theory,  which  mitigates  the 
difficulties  of  the  previous  theories,  are  as  follows:  There  must  be 
a  mechanism  to  generate  a  sufficient  number  of  starting  electrons  to 
initiate  the  avalanche.  These  few  starting  electrons  drift  upward  in 
energy  in  the  conduction  band  through  interaction  with  the  laser  field, 
retarded  by  energy  loss  to  phonons.  Each  electron  undergoes  a  random 
walk  of  progressively  increasing  kinetic  energy  until  it  attains  a 
threshold  energy  £j.  The  electron  then  drops  to  the  bottom  of  the 
conduction  band.  If  an  exciton  has  been  created,  its  rapid  photo¬ 
ionization  injects  a  second  electron  into  the  bottom  of  the  conduction 

band.  The  process  is  repeated  until  a  sufficient  electron  density 
18  -3 

(~10  cm  )  has  been  created  to  damage  the  crystal,  probably  by  ex¬ 
cessive  Joule  heating. 

This  description  of  electron-avalanche  breakdown  is  indeed  also 
common  to  most  previous  theories.  However,  a  cardinal  reason  for  the 
good  agreement  of  the  present  theory  with  experimental  results,  con¬ 
trasted  with  the  poor  agreement  of  previous  similar  theories,  is  the 
use  of  more  realistic  electron-phonon  scattering  rates.  Umklapp  pro¬ 
cesses,  which  were  neglected  in  previous  treatments'6’^7  are  included 
in  the  electron-phonon  relaxation  frequencies.  In  other  words,  the 
phonons  with  which  the  electrons  interact  are  not  restricted  to  those 
in  the  first  Brillouin  zone,  as  in  previous  theories.  Both  acoustical 
and  optical  phonons  are  included  in  the  interactions  with  electrons. 

Our  theory  retains  the  difference  between  the  electron  relaxation  fre¬ 
quencies  (or  scattering  rates)  yk  at  which  the  electron  loses  its 
forward  component  of  momentum  and  yL  at  which  the  electron  loses 
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energy.  The  resulting  magnitude  and  energy  dependences  of  the  electron- 
phonon  scattering  rates  are  essential  differences  between  our  model  and 
those  of  previous  investigations. 

Another  unique  aspect  of  the  present  theory  is  the  inclusion  of 
the  temperature  dependence  of  the  lattice  constant  and  phonon  fre¬ 
quencies.  Without  that  dependence,  the  agreement  with  the  measured 
temperature  dependence  of  the  breakdown  irradiance--especially  at 
10.6  ym--is  degraded.  With  the  inclusion  of  the  new  electron  scatter¬ 
ing  rates  and  the  temperature  dependence  of  the  parameters,  the  common 
transport-equation  approach  can  explain  available  data  without  intro¬ 
ducing  new  mechanisms.  In  particular,  the  predictions  of  the  theory 

13  17 

agree  well  with  the  new  Lebedev  data.  * 

The  material -to-material  variations  of  the  breakdown  electric 
field  Eg  in  the  present  theory  result  from  the  values  of  the  character¬ 
istic  phonon  frequency  which  is  set  equal  to  the  Debye  energy 

kgeD,  and  the  electron-cation  scattering  cross  section  Qa>  which  is 
scaled  as  the  square  of  the  ionic  radius. 

The  electrons  that  are  initially  present  or  generated  (those  that 
start  the  avalanche)  will  be  called  the  starting  electrons.  Most 
previous  theories  of  electron-avalanche  breakdown  in  solids  at  laser 
frequencies  have  either  implicitly  or  explicitly  assumed  that  starting 
electrons  are  present  when  the  laser  is  turned  on.  Sparks17  showed 
that  experimental  data  contradict  that  assumption.  The  starting  elec¬ 
trons  are  not  initially  present;  they  must  be  generated--by  an  external 
source  or  by  the  laser  field  itself,  for  example.  In  microwave  gas- 
breakdown  experiments,  it  is  well  known  that  the  starting  electrons 


0 


-5- 

must  be  generated.  A  radioactive  source  near  the  microwave  cavity  is 
commonly  used  to  provide  the  starting  electrons. 

Since  both  processes--generation  and  acceleration--are  required 
for  breakdown,  either  the  generation  of  the  starting  electrons  or  the 
avalanche  process  itself  can  determine  the  breakdown  field.  Processes 
by  which  the  laser  field  can  generate  the  starting  electrons  were  dis¬ 
cussed  by  Sparks.17  Only  the  avalanche  process  is  considered  here. 

The  estimates  of  Sparks17  on  electron  generation— as  well  as  the  good 
agreement  between  experiment  and  the  present  theory,  in  which  the 
breakdown  is  dominated  by  the  avalanche  process— partially  justify  not 
conducting  a  detailed  investigation  of  starting-electron  generation. 

In  fact,  the  importance  of  generating  the  starting  electrons  will  re¬ 
main  uncertain  until  verified  experimentally. 

If  the  electrical  breakdown  field  Eg  were  dominated  by  starting 
electrons  either  initially  present  or  generated,  rather  than  by  the 
avalanche  process  itself,  the  experimental  value  of  Eg  would  generally 
be  greater  than  that  obtained  when  starting  electrons  are  provided. 

If  dominated  by  electrons  initially  present.  Eg  would  vary  from  shot 
to  shot,  depending  on  the  probability  of  a  starting  electron  being 
present.  Such  variations  are  well  known  in  microwave  breakdown  of 
gases,  in  which  Eg  can  vary  by  well  over  an  order  of  magnitude. 

A  major  difficulty  in  obtaining  agreement  between  theory  and 
experiment  is  determining  whether  the  experimental  values  of  Eg  are 
intrinsic.  It  is  common  practice  to  select  only  the  greatest  mea¬ 
sured  values  of  Eg.  However,  those  retained  still  may  not  be  the 
intrinsic  values.  Despite  the  inherent  difficulty  of  determining 


when  an  intrinsic  value  has  been  attained,  it  is  likely  that  the 
Lebedev  values  at  1.06  and  10.6  pm  are  intrinsic  because  they  agree 
so  well  with  the  present  theory. 

In  using  a  differential  form  of  the  transport  equation,  we  assume 
that  the  photon  energy  fa  is  sufficiently  small  to  be  treated  by  dif¬ 
ferentials  rather  than  differences.  At  high  frequencies  (wavelength 
X  <  1  pm,  very  roughly),  the  differential  approximation  gives  rise  to 

nonnegligible  errors.  This  and  other  high-frequency  effects  will  be 

19 

considered  in  a  subsequent  publication. 

Sparks^  showed  that  the  process  in  which  an  electron  absorbs  or 

emits  both  a  photon  and  a  phonon  is  sometimes  important  in  electron- 

avalanche  breakdown.  This  was  called  the  Holstein  process  because 
20 

Holstein  showed  it  to  be  the  dominant  mechanism  determining  the 
value  of  the  conductivity  in  the  limit  ojt  »  l--provided  that  fa  is 
sufficiently  small  to  permit  photon  emission  as  well  as  absorption. 

The  Holstein  process  is  tacitly  included  here  because  the  energy  gain 
of  the  electron  is  treated  in  terms  of  the  conductivity  (which  is  con¬ 
trolled  by  the  Holstein  process  in  the  limit  urr  »  1). 

Table  I  lists  processes  not  included  in  the  present  theory  be¬ 
cause  they  are  believed  to  be  negligible  in  careful  experiments  for 
h  1.06  pm,  along  with  the  expected  effect  on  the  value  of  Eg. 

The  major  sections  of  this  paper  are  as  follows: 

I.  Introduction 

II.  Average-Electron  Model 

III.  Diffusion  Model 

IV.  Electron  Relaxation  Frequencies 


m 
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V.  Numerical  Results 

VI.  Interpretations  and  Calculations  Using  Simple  Models 
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II.  AVERAGE-ELECTRON  MODEL 


In  order  to  understand  the  diffusion  model,  which  is  the  basis 
of  the  present  investigation,  first  consider  the  average-electron 
model.  Section  I  discussed  the  general  features  of  most  current  model 
the  acceleration  of  a  conduction  electron  by  the  laser  field,  the  loss 
of  energy  from  the  electron  to  the  phonons,  the  generation  of  a  second 
conduction  electron  accompanied  by  a  loss  of  kinetic  energy  of  the 
first  electron,  and  the  repetition  of  the  process  until  the  electron 
concentration  is  sufficient  to  damage  the  crystal. 

In  the  average-electron  model,  which  affords  the  simplest  mathe¬ 
matical  treatment  of  the  general  model,  the  average  energy-loss  rate 
to  the  lattice  is  (d£/dt)E  and  the  average  energy-gain  rate  from  the 
electric  field  is  (d£/dt)E-  Breakdown  occurs  when  the  gain  exceeds 
the  loss: 

(d£/dt)E  >  (d£/dt)L  (2.1) 

for  all  &  <  where  is  the  energy  at  which  the  electron  generates 
a  second  conduction  electron  by  excitation  across  the  electronic 
energy  gap. 

The  rate  at  which  the  field  £  adds  energy  to  the  electron  is  ob¬ 
tained  from  the  wave-vector  equation 

fl{ dkz/dt  +  Ykkz)  =  eEp|<  exp  (iwt)  ,  (2.2) 

where  E,  is  along  the  z  axis,  kz  is  the  z  component  of  the  electron 
wave  vector  k,,  y^  =  l/xk  is  the  transport  (momentum  loss)  relaxation 


■■ 
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frequency,  and  u  is  the  laser  frequency.  The  steady-state  solution 
is 

kz  ■  ffll  +\T|<)  6pt  «»p  0-0  ■  <2'3> 

With  current  ^  =  ey,  where  the  2  component  v2  of  the  velocity  y,  is 
equal  to  and  m*  is  the  electron  effective  mass,  the  rate  at 

which  the  field  puts  energy  into  the  electron  is 

(d£/dt)r  =  i  Re  J  •  E*  =  l  Re  fie k,E*/m*  , 

t  c  **  +*  c  Z 

where  Re  denotes  the  real  part  and  Ej*  is  the  complex  conjugate  of  E. 

2  12 

Substituting  Eq.  (2.3)  into  this  equation,  using  E  =  ^  Epk  (where 
E  is  the  root-mean-square  value  of  the  field),  and  taking  the  average 
over  the  electron  distribution,  we  obtain 

e2r 

(d<S/dt)E  =  1  oE2  ,  a  *  - K-j-  ,  (2.4) 

J  m*(l  +  w  t£) 

where  a  is  the  electrical  conductivity.  Equation  (2.4)  is  the  well- 
known  conductivity  result  for  a  single  electron.  Physically,  the 
energy  gain  in  Eq.  (2.4)  is  small  when  the  electron  momentum  is  changed 
rapidly  by  interactions  with  the  phonons.  That  is,  in  this  dc  limit 
of  large  (or,  more  precisely,  small  iot^),  the  phonon  collisions 
inhibit  the  acceleration  of  the  electron  by  changing  its  momentum  to 
the  opposite  direction  from  that  imparted  by  the  electric  field. 

In  the  opposite  limit  of  slow  momentum  change,  the  energy  gain 
in  Eq.  (2.4)  is  again  small.  For  this  case  of  »  1,  the  electron 
is  first  accelerated,  then  deaccelerated  when  the  direction  of  the 
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electric  field  reverses.  The  process  occurs  many  times  before  the 
electron  undergoes  a  momentum- chang inq  collision  with  a  phonon.  The 
energy  gained  during  one  helf-cycle  is  lost  during  the  following  half¬ 
cycle.  The  maximum  energy  gain,  which  is  a  function  of  the  momentum 
relaxation  time  t^,  occurs  at  uvr^  =  1;  then,  on  the  average,  the  elec¬ 
tron  direction  is  reversed  every  time  the  electric  field  reverses  direction. 

The  rate  at  which  the  electron  loses  energy  to  the  lattice  is 

(dS/dt)L  =  &)phYL  ,  (2.5) 


where  uiph  is  the  average  phonon  frequency  and  the  energy-loss  relaxa¬ 
tion  frequency  yl  is  different  from  yk  in  Eq.  (2.4),  as  discussed 
below.  Equating  Eqs.  (2.4)  and  (2.5)  gives  the  threshold  value  EVQ 
above  which  the  electrons  gain  energy  from  the  field 


'VO 


-ewrM-C 


(2.6) 


The  resulting  breakdown  criterion  that  the  electric  field  E  must  be 
greater  than  EVQ  at  all  electron  energies— that  is. 


EB  =  EV0mx 


(2.7) 


where  Eygmx  is  the  maximum  value  of  E^q,  and  Eg  is  the  value  of  E  at 
breakdown— is  a  fair  zeroth-order  approximation,  provided  the  energy- 
dependent  relaxation  frequencies  y^  and  yk  are  interpreted  correctly. 
However,  even  with  the  proper  interpretation  of  the  relaxation  fre¬ 
quencies,  the  average-electron  result  in  Eqs.  (2.6)  and  (2.7)  is  only 


a  rough  approximation.  As  discussed  below,  the  inadequacy  of  the 
average-electron  result  stems  both  from  the  difficulty  of  treating 


the  energy  dependence  y^  and  y^  properly,  and  from  the  neglect  of 
electron  diffusion  in  energy  space.  The  problems  are  overcome  by 
using  the  transport-equation  approach  discussed  in  the  following 
section. 

Before  leaving  the  average-electron  model,  consider  the  multipli¬ 
cation  rate  3.  As  a  result  of  the  repeated  doubling  of  the  number 
of  conduction  electrons,  their  number  Ng  increases  exponentially  as 
Ne  =  N0  exp  (8*)*  The  net  energy  gain,  from  Eqs.  (2.4)  and  (2.5),  is 

_  e^ 

<*dt  -  i  oE2  -  =  1  oE20  .  (2.8) 

Lvo 

o 

Assuming  that  a  and  EyQ  are  independent  of  t  (that  is,  yk  and  y^ 
independent  of  £) ,  integrating  Eq.  (2.8),  and  using  the  breakdown 
criterion  of  setting  the  pulse  duration  t  =  tp  equal  to  51/8  as  dis¬ 
cussed  in  Sec.  V,  we  obtain 


8  = 


51  E2  -  E2q 


for  E  >  E 


’I 


'VO 


=  0  , 


for  E  <  E 


VO  ’ 

2 

VO  • 


(2.9) 


The  value  of  8  goes  to  zero  at  E  =  Eyg.  For  a  high  electric  field 
E2  »  E2q,  Eq.  (2.9)  gives 

__1  0  O  O 

(2.10) 


8  =  17£j]aE2  ,  for  E2  »  E2Q  . 


— 
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III.  DIFFUSION  MODEL 


The  central  approach  of  the  present  investigation  is  to  solve 
the  appropriate  transport  equation  by  an  eigenvalue  method.  The  trans¬ 
port  equation,  detailed  by  Holway  and  FradinJ®  can  be  cast  as  a  dif¬ 
fusion  equation  in  energy  space: 

3n(2,  t)/3t  +  3J/3&  =  0  ,  (3.1) 

where  n(&,  t)  is  the  number  of  electrons  with  energy  between  t  and 
£+  d£,  and  the  energy-space  current  (the  net  number  of  electrons 
whose  energy  increases  from  a  value  less  than  L  to  a  value  greater 
than  &  per  unit  time)  is 

J(£,  t)  =  V(£)n(£,  t)  -  D(S)9n(S,  t)/dZ  .  (3.2) 

Here  V(&)  and  D(£)  are  the  effective  velocity  and  diffusion  coefficient 

in  energy  space,  respectively.  Both  depend  on  energy,  temperature,  and 

✓ 

laser  frequency.  For  fields  comparable  to  the  breakdown  field,  ignoring 
small  field-independent  contributions  to  D,  we  obtain 

V(£)  =  joE2  -  4)phYL  ,  D(<e)  =  |  &E2  ,  (3.3) 

where  E  is  the  root-mean-square  electric  field  and  o  is  defined  in 
Eq.  (2.4).  The  boundary  conditions  are 


n(£j,  t)  =  0 


(3.4) 
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1 


corresponding  to  no  electrons  accelerated  past  the  exciton  threshold, 
and 

J(0,  t)  =  t)  ,  (3.5) 

corresponding  to  the  "flux  doubling"  feature  of  exciton  creation 
followed  by  photoionization.  Converting  the  transport  equation  to  a 
diffusion  equation  [Eq.  (3.1)]  is  extremely  useful  because  well-known, 
intuitively  understood  diffusion  results  can  be  brought  to  bear  on  the 
electron-avalanche  breakdown  problem,  as  illustrated  below. 

The  transport  equation  resulting  in  Eqs.  (3.1)  through  (3.5)  has 

20 

been  derived  in  detail  by  Holstein  --starting  with  the  Boltzmann 
equation,  neglecting  the  spatial  dependence  of  the  electron  distribu¬ 
tion,  and  averaging  over  momentum  directions  to  reduce  the  momentum 
dependence  to  an  energy  dependence.  In  starting  with  the  Boltzmann 
differential  equation,  it  is  tacitly  assumed  that  the  photon  energy 
hn  and  the  phonon  energies  h are  sufficiently  small  that  energy 
differences  can  be  replaced  by  differentials.  The  transport-equation 
results  in  Eqs.  (3.1)  through  (3.5)  can  also  be  obtained  from  a 
straightforward  summation  of  transition  probabilities. 

The  partial -differential  equation  (3.1)  is  converted  to  an 
ordinary-differential  eigenvalue  equation  by  using  the  approximation 

n(2,  t)  =  n(£)  exp  (3t)  .  (3.6) 

Substituting  Eq.  (3.6)  into  Eq.  (3.1)  gives  the  eigenvalue  equation 


dJ(£)/d£  =  3n(&)  . 


(3.7) 
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The  eigenvalue  6  is  the  multiplication  rate,  and  the  eigenfunction 
n(&)  is  the  electron  density  in  energy  space.  When  the  laser  is  first 
turned  on,  the  initial  electron  distribution  n(£,  0)  is  not  of  the  form 
assumed  in  Eq.  (3.6).  However,  the  initial  distribution  quickly  de¬ 
velops  into  the  eigenfunction  distribution  n(£)  in  a  time  of  the  order 
of  the  diffusion  time  constant 

Tdif  =  l\,U  *  (3.8) 

In  fact,  for  the  exactly  solvable  example  of  V  =  0  and  D  constant,  the 
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approach  time  to  n(£)  is  a  factor  of  ten  shorter  than  t^.  Since 
Tdif  1S  ln  t^ie  subPicosecond  region  in  cases  considered  to  date,  the 
form  n(&)  is  attained  over  essentially  the  full  pulse  duration,  and 
the  approximate  eigenvalue  solution  is  quite  accurate. 

The  accuracy  of  the  eigenvalue-equation  solution  has  also  been 
demonstrated  explicitly  by  expanding  a  representative  initial  electron 
distribution  n(£,  0)  in  the  eigenfunctions.  Then  the  projection  of 
n(£,  0)  on  the  eigenfunction  with  positive  real  B  (as  discussed  below) 
grows  according  to  Eq.  (3.6)  for  all  t  >  0,  and  the  projections  on 
the  other  eigenvectors  decay  rapidly.  Dual-space  eigenfunctions  must 
be  used  in  the  scalar  products  with  which  the  expansion  coefficients 
are  computed,  because  the  system  operator  is  non-Hermitian. 

An  eigenvalue  B  with  a  positive  real  part  corresponds  to  an  elec¬ 
tron  concentration  that  increases  exponentially  in  time  [according  to 
Eq.  (3.6)],  as  expected  for  the  avalanche  process.  The  eigenvalue 
problem  defined  by  the  differential  equation  (3.7),  with  boundary 
conditions  in  Eqs.  (3.4)  and  (3.5),  is  non-Hermitian  and  nonpositive 


definite;  thus  there  is  no  assurance  that  either  real  eigenvalues  or 
positive  eigenvalues  exist,  though  the  eigenvalues  is  must  occur  in 
complex-conjugate  pairs.  Insight  is  gained  from  three  simple  cases 
with  exact  solutions  to  the  eigenvalue  problem: 

1.  V(£)  =  -Vq  =  constant  and  D(<£)  =  Dq  =  constant. 

2.  V(£)  =  -Vq  =  constant,  with  Vq  >  0,  and  D(£)  =  D^£. 

3.  V(&)  =  Vq  -  V,£  and  D(£)  =  2 VQ£. 
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These  cases  are  discussed  elsewhere.  In  all  three,  there  is  one  and 
only  one  eigenvalue  with  a  positive  real  part  (when  E  is  sufficiently 
large),  and  the  imaginary  part  of  that  eigenvalue  is  zero.  All  others 
have  Re(B)  <  0,  which  corresponds  to  an  exponentially  decreasing  elec¬ 
tron  concentration.  In  our  numerical  work,  we  have  always  found  one 
and  only  one  eigenvalue  on  the  positive  real  axis.  We  suspect  that 
this  is  a  general  feature  of  the  present  problem,  but  we  have  no  proof. 

The  boundary  conditions  in  Eqs.  (3.4)  and  (3.5)  are  unrealistic 
in  one  regard.  After  creation  and  subsequent  photoionization  of  the 
exciton,  both  electrons  cannot  appear  together  at  precisely  zero  energy, 
but  are  in  fact  injected  into  the  band  energy  over  an  energy  interval 
with  finite  width  controlled  by  both  the  exciton  linewidth  (for  sodium 
chloride  at  room  temperature,  a  few  tenths  of  a  volt)  and  the  energy 
dependence  of  the  photoionization  cross  section.  In  our  numerical 
work,  we  have  explored  the  sensitivity  of  the  solution  to  this  boundary 
condition  by  modifying  Eq.  (3.6)  with  a  source  term  that  distributes  a 
total  of  two  electrons  over  a  finite  energy  interval  near  the  bottom 
of  the  conduction  band,  for  each  electron  that  strikes  the  threshold 
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for  exciton  creation  at  £j.  The  solution  is  affected  only  slightly 
by  that  procedure,  even  when  the  "injection  width"  is  several  tenths 
of  an  electron  volt. 


IV.  ELECTRON  RELAXATION  FREQUENCIES 


The  transport  relaxation  frequency  yk  and  the  energy-loss  param¬ 
eter  /to phYL  are  obtained  from  the  expressions 


Yk  =  I  I  E  (1  -  k  •  k')pq±(k,  k‘) 


a.  K. 


(4.1) 


and 


ShYL  =  N(£)'1 


nn 

±qkk' 

A*  #**  rw 


(±\)pq±(il»  -  &k)  . 

•  “  iv 


(4.2) 


where  the  caret  (-)  denotes  a  unit  vector,  pq+(k.»  k,1 )  is  the  proba¬ 
bility  that  an  electron  is  scattered  from  k,  to  k,1  by  absorption  or 
emission  (±)  of  a  phonon  of  wave  vector  q  and  frequency  u>  ,  6  is  the 

q 

Dirac  6  function,  and  N(£)  is  the  density  of  states  in  the  conduction 
band.  The  summations  over  include  summations  over  phonon  branches. 
The  expression  in  Eq.  (4.2)  provides  the  value  of  the  rate  of  energy 
loss  ^p^Yl  directly,  rather  than  y^.  A  value  for  yL  itself  can  be 
obtained  by  setting  the  average  phonon  energy  equal  to  the  Debye 
energy 

^ph  =  kB6D  *  (4-2a) 

where  8Q  is  the  Debye  temperature  and  kg  is  the  Boltzmann  constant. 

First  consider  the  contribution  of  the  longitudinal-optical 
phonons  to  y^.  In  polar  materials  such  as  the  alkali  halides,  the 
electrons  couple  strongly  to  long-wavelength  longitudinal -optical 
phonons  as  a  result  of  the  macroscopic  electric  field  they  generate. 
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Within  the  framework  of  a  model  in  which  the  electron  moves  in  a  para¬ 
bolic  band  with  effective  mass  m*,  it  is  straightforward  to  derive 
the  contribution  to  y^  from  Frohlich  scattering,  as  the  scattering 
from  the  longitudinal-optical  phonons  is  called.  When  the  electron 
energy  &  is  large  compared  with  the  longitudinal-optical  phonon  fre¬ 
quency  the  result  of  evaluating  the  summations  in  Eq.  (4.1)  is^7 


2  *1/2 
e  m*  u), 


YkL0 


10 


Jim 


i*V2 


(ti) 


(1  +  2nLQ)  , 


(4.3) 


where  nLQ  =  [exp  (/to^Q/kgT)  -  1]“"'  is  the  number  of  thermal  longi¬ 
tudinal-optical  phonons  present  in  thermal  equilibrium  at  temperature 
T,  and  €q  and  6^  ere  the  static  and  high-frequency  dielectric  constants, 
respectively. 

The  corresponding  result  for  ^phYL  from  Frohlich  scattering  is 


%h\  - 


V2  2  2 

(2m*)  e  WL0  /  1  1 


£ 


1/2 


te- *)-(£) 


From  the  <T172  dependence  displayed  in  Eqs.  (4.3)  and  (4.4),  we  see 
that  Frohlich  scattering  is  most  important  near  the  bottom  of  the 
conduction  band.  The  behavior  displayed  in  Eqs.  (4.3)  and  (4.4) 

is  valid  only  when  £  »  /&>Lq.  The  full,  complex  expression  is  approx¬ 
imated  by  rounding  off  the  singularity  by  replacing  £  with  2/a^g  for 
electron  energies  <  2 This  approximation  has  little  influence 
on  the  results. 

Next  consider  the  contribution  to  y^  from  nonpolar  phonons--that 
is,  from  acoustical  and  transverse-optical  phonons.  The  probability 
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that  an  electron  is  scattered  from  to  k,'  by  emission  or  absorption 
of  a  phonon  of  wave  vector  £  and  frequency  is  given  by 

Pq+(ji*  L')  =  X  Mkk *  |2  (mLt)  ~  ~  + 


x  (nq  +  1/2  +  l/2)6(«Sk,  -  \  +  /<oq) 


(4.5) 


where  M  is  the  effective  ionic  mass,  N  is  the  number  of  unit  cells  per 
unit  volume,  A  is  the  Kronecker  delta,  6  is  the  Dirac  <5  function,  and 
Mkk,  is  the  electron-phonon  matrix  element  defined  below.  For  acous¬ 
tical  phonons  of  long  wavelength,  M  is  the  sum  of  the  masses  in  the 
unit  cell;  for  transverse-optical  phonons  of  long  wavelength,  M  is  the 
reduced  mass  of  the  unit  cell;  and  for  zone  boundary  phonons,  M  is 
a  weighted  average  of  the  mass  on  each  sublattice.  Since  the  extended 
zone  scheme  is  used,  ^  is  not  restricted  to  the  first  Brillouin  zone 
and  the  phonon  frequency  is  a  periodic  function  of  the  wave  vector, 
with  the  periodicity  of  the  Brillouin  zone. 

Within  the  framework  of  the  simplest  possible  theory--the  free- 

electron-diffraction  theory  with  the  Born  approximation--the  electron- 
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phonon  matrix  element  is 


Mkk'  ■  iVa(K)(K  •  eQ)  , 


(4.6) 


where  K  =  k,  -  k/ ,  e  is  the  unit  phonon  eigenvector,  and 


=  na/^P  (“  ™  • 


(4.7) 
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where  na  is  the  atomic  density  and  l)  (r)  is  the  potential  in  a  unit 
cell  of  the  crystal.  The  validity  of  Eqs.  (4.6)  and  (4.7)  can  be  ex¬ 
tended  beyond  the  Born-approximation  limit  by  eliminating  |v  (K)l^ 
through  use  of  the  expression 


for  the  cross  section  per  unit  solid  angle  a  (K). 

a 

For  sodium  chloride,  examination  of  the  magnitude  of  the  cross 
sections  for  scattering  from  Na+  and  Cl"  ions  shows  that,  in  the  range 
of  energies  of  interest  here,  the  cross  section  for  scattering  from 
the  large,  polarizable  chlorine  is  far  greater  than  that  of  the  sodium 
ion.  The  corresponding  result  is  true  for  all  alkali  halides.  Thus, 
o  ( K )  is  well  approximated  by  the  cross  section  for  scattering  by  the 

a 

negative  ion.  In  the  numerical  calculations,  the  value  of  cr  (K)  is 

d 

taken  as  the  free-ion  cross  section.  In  effect,  that  procedure  extends 
the  validity  beyond  the  Born  approximation  because  the  full  t-matrix 
treatment  for  scattering  by  the  negative  ion  is  used  in  the  evaluation 
of  a (K). 

d 

For  values  of  q  outside  the  Brillouin  zone,  the  scattering  events 
described  by  Eq.  (4.5)  are  Umklapp  processes.  A  proper  description  of 
Umklapp  scattering  requires  a  more  detailed  calculation  of  Mkk-  Even 
though  the  full  t-matrix  for  scattering  by  the  halide  ion  is  used  to 
evaluate  a,(K),  the  calculation  in  effect  treats  the  electron  wave 

a 

function  using  a  one-plane-wave  approximation.  This  leads  to  a  non¬ 
physical  divergence  in  Pq+(k,,  k^‘ )  whenever  the  wave-vector  transfer 
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equals  a  reciprocal  lattice  vector  £  because--as  £  approaches  G  for 

the  three  acoustical  branches  of  the  phonon  spectrum--^  approaches 

zero  while  M^,  remains  finite.  A  proper  treatment  of  the  electron- 
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phonon  matrix  element,  which  is  equivalent  to  at  least  a  two-plane- 
wave  approximation  of  the  matrix  element,  shows  that  the  electron- 
phonon  matrix  element  vanishes  as  £’  approaches  £  +  £,  and  that  the 
cross  section  varies  smoothly  near  these  special  values  of  wave-vector 
transfer. 

When  £  lies  outside  the  first  Brillouin  zone,  the  phonon  fre¬ 
quencies  are  approximated  with  suitable  weighted  averages  of  zone¬ 
boundary  phonon  frequencies,  and  an  average  is  taken  over  the  direc¬ 
tions  of  the  phonon  eigenvector  e^.  That  approximation  properly  in- 

r\ 

eludes  the  effect  of  the  factor  (K  •  e  )*  on  the  energy  variation  of 
the  cross  section  as  well  as  the  correct  volume  of  phase  space  avail¬ 
able  to  the  electron  in  its  final  state.  The  approximated  value  of 
the  cross  section  then  has  the  correct  overall  energy  dependence,  with 
a  reliable  semiquantitative  estimate  of  its  magnitude  at  large  electron 
energies  where  Umklapp  scattering  dominates. 

The  limit  of  small  electron  wave  vectors  is  k  ~  kgZ/2  and 
£  «  z/2,  where  kgZ  is  the  wave  vector  of  an  electron  at  the  Brillouin 

zone  boundary.  If  we  use  that  limit  in  evaluating  the  summations  in 
Eq.  (4.1)— with  only  the  summations  over  acoustical  phonons  included, 
because  the  available  phase  space  with  q  <  2k  is  small—we  obtain 


Y  k(£)  =  A ^-5 

8w/ 


2k 


C  4  ftCAq 

J  dqq  coth  2y  * 


for  &<  j  &BZ 


(4.9) 
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where  ag  is  the  lattice  constant,  is  the  total  mass  of  the  unit 
cell,  c^  is  the  sound  velocity,  k  is  the  wave  vector  of  the  electron, 
and  Q(£)  =  4ira  (6)  is  the  total  cross  section  integrated  over  all 

a 

solid  angles.  In  the  limit  of  small  k,  Eq.  (4.9)  gives 


1 


V 


/9m*vl/ 2  3.,  2 

(2m*)  a0Mrc£ 


&1/2Q(dO 


(4.10) 


for  a  parabolic  band  with  effective  mass  m*. 

In  the  opposite  limit  k  »  k^/2,  the  electron  effective  mass 
is  approximated  with  the  free-electron  mass  mg  and  the  corresponding 
result  is 


?  ?  /  \ 3/2 

Yk(£)  =  ~ -6-9-(-41  ( \ 
a0M>m0  '^BZ  / 


(1  +  2nu)  (1  +  2nTfl) 


0), 


+  2 


TA' 


LA 


to 


TA 


(4.11) 


where  taLA  and  tjyA  are  the  frequencies  of  longitudinal  and  transverse 
phonons  at  the  Brillouin  zone  boundary,  respectively,  and  the  effec¬ 
tive  mass  M  is  chosen  as  M>--the  heaviest  constituent  mass  in  the  unit 
cell . 

In  the  numerical  calculations,  we  use  a  straightforward  extrap¬ 
olation  between  the  low-energy-dominated  relaxation  frequency  in  Eq. 
(4.9)  and  the  high-energy,  Umklapp-dominated  relaxation  frequency  in 
Eq.  (4.11).  The  high-energy  behavior  of  y^(£)  displayed  in  Eq.  (4.11) 

is  a  crucial  feature  of  our  analysis.  Both  Holway  and  Fradin^  and 
12  13 

the  Soviet  analysts  ’  neglected  Umklapp  scattering,  thereby  missing 
the  dominant  contribution  to  the  relaxation  frequencies  for  electron 
energies  high  in  the  conduction  band. 
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Holway  and  Fradin  approximated  the  acoustical  phonon  scattering 
with  an  expression  similar  to  Eq.  (4.9),  but  replacing  the  upper  limit 
on  the  integral  with  a  constant  for  the  case  of  i  Us’nf)  that 

constant,  high-energy  upper  limit  is  equivalent  to  including  only 
normal  processes  in  the  expression  for  vk(<5).  In  that  case,  the 
scattering  rate  decreases  with  energy  at  high  energies--in  qualitative 
contrast  to  Eq.  (4.11) 

The  Soviet  work  employed  an  expression  equivalent  to  Eq.  (4.10), 
but  with  Frohlich  scattering  ignored.  The  resulting  form  for  Yk(£) 
was  extrapolated  to  energies  beyond  £gZ,  thereby  underestimating  the 
strength  of  the  acoustical  phonon  scattering  at  high  energies. 

Finally,  consider  the  contribution  to  ft* phY|_  from  nonpolar  scat¬ 
tering.  In  calculating  that  contribution,  the  phonon  frequencies  that 
appear  in  the  Dirac  6  function  in  Eq.  (4.5)  are  neglected,  and  approx¬ 
imations  similar  to  those  above  are  used.  For  scattering  from  a  single 
phonon  branch,  this  yields 

n  v(£)/t2  r  „  9 

^ph^L  =  2M  J  k  )^a(K)(K  •  e^)  ,  (4.12) 

where  v(6)  is  the  velocity  of  an  electron  of  energy  o.  In  the  low- 
energy  limit,  where  only  long-wavelength  acoustical  phonons  contribute, 
evaluating  the  integral  in  Eq.  (4.12)  gives 

/*>phyL  =  2nav(S)Q(£)  .  (4.13) 

For  t  >  £gZ/2,  m*  is  replaced  with  the  free-electron  mass  and  Mj. 
is  replaced  with  M>  .  The  numerical  calculation  of  the  Frohlich 
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contributions  to  Yk  and  /<WphY|_  is  straightforward.  Only  the  value 
of  m*  near  the  bottom  of  the  conduction  band  is  needed.  The  usual 
approximation  m*  s-  itIq/2  is  used,  where  is  the  free-electron  mass. 

To  calculate  the  nonpolar  phonon  contributions  to  y^  and  y^*  the 
cross  section  Q  (<S)  for  scattering  of  an  electron  by  the  halide  ion 

d 

is  replaced  with  its  value  Qa(£j)  at  the  exciton  energy.  Then  Qg(£j) 

is  estimated  from  molecular  scattering  data.  For  example,  the  cross 

section  for  electron  scattering  from  the  molecule  HC1  is  dominated  by 

scattering  from  the  Cl’  ion.  Thus,  we  approximate  Qa(oj)  with  the  HC1 

2 

cross  section,  which  has  the  value  Q(£j)  a*  0.35  nm  at  8  eV. 

2 

For  the  other  alkali  halides,  we  use  the  scaling  Q(£j)  =  ( rh/rNaCl ^ 

2 

0.35  nm  ,  where  rh  is  the  ionic  radius  of  the  halide  ion.  The 
estimated  values  of  Q(«Sj)  for  sodium  chloride  and  the  other  alkali 
halides  are  some  of  the  major  uncertainties  in  obtaining  quantitative 
values  for  the  acoustical  phonon  contributions  to  Yk  and  yl*  The 
phonon  frequencies  are  approximated  using  =  c-j-k^  and  u>|^  = 
c£kBZ’  for  q  ^  kBZ’  wliere  kBZ  1S  tlie  diameter  of  a  spherical  Brillouin 
zone  chosen  to  have  volume  equal  to  the  real  Brillouin  zone. 

The  temperature  dependence  of  the  relaxation  frequencies  is  de¬ 
termined  both  by  the  explicit  temperature  dependence  of  the  Bose- 
Einstein  phonon-occupation  numbers  and  by  the  temperature  dependence 
of  the  lattice  constant  and  the  phonon  frequencies.  Figures  1  and  2 
show  the  magnitude  and  energy  variations  of  Yk  and  Y|_  for  sodium 
chloride  at  room  temperature.  The  Frohlich  contribution  and  the 
acoustical  phonon  contribution  are  shown  separately  to  illustrate  their 
relative  importance.  The  Frohlich  contribution  is  important  only  at 
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low  energies,  within  one  to  two  volts  of  the  band  edge.  Fioure  3 
shows  the  full  curves  of  and  Y|_  together. 

Values  of  the  parameters  used  in  the  numerical  calculations  are 
listed  in  Table  II. 


^«si 
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V.  NUMERICAL  RESULTS 


The  electron  multiplication  rate  6  was  calculated-using 
Hamming’s  predictor-corrector  method— by  a  careful  numerical  inte¬ 
gration  of  the  eigenvalue  equation  in  Eq.  (3.7),  with  boundary  con¬ 
ditions  defined  by  Eqs.  (3.4)  and  (3.5).  If  we  consider  the  behavior 
of  3  as  a  function  of  electric  field  and  temperature,  then  the  break¬ 
down  fields  predicted  by  the  model  will  be  obtained  from  the  3(E) 
curves.  The  dependence  of  3  on  the  electric  field  E  is  shown  in  Fig. 

4  for  sodium  chloride  at  1.06  ym  and  four  temperatures.  Similar  curves 

for  10.6  ym  are  shown  in  Fig.  5.  The  dramatic  increase  in  3  with  in- 

2  2 

creasing  E  at  low  fields,  E  «  Eygmx,  is  explained  in  Sec.  VI  and  can 

be  seen  from  the  average-electron  model,  for  which  3  goes  exactly 

to  zero,  as  in  Eq.  (2.9). 

2  2 

At  high  fields,  E  »  E VOmx,  the  multiplication  frequency  has  the 
general  behavior 

3  ~  E2  ,  for  E2  »  E2q  ,  (5.1) 

as  in  Eq.  (2.10)  and  as  discussed  in  Sec.  VI.  The  curves  in  Figs.  4 

and  5  do  not  extend  to  sufficiently  great  electric  fields  to  fully 

2  2 
show  the  E  dependence,  but  we  have  verified  that  the  quite  general  E 

dependence  is  attained  at  high  fields.  The  temperature  variation  of 

3  displayed  in  Figs.  4  and  5  can  be  understood  qualitatively,  as  also 

discussed  in  Sec.  VI. 

The  breakdown  field  is  easily  obtained  by  calculating  the  multi¬ 
plication  frequency  3  as  a  function  of  electric  field  once  a  criterion 
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has  been  established  for  selecting  the  value  of  6.  We  use  the  criterion 
8tg^^  =  17,  where  tg^  is  the  effective  pulse  duration.  For  Gaussian 
pulses,  teff  is  set  equal  to  one-third  the  total  pulse  duration  be¬ 
cause  breakdown  occurs  principally  near  the  center  of  the  pulse.  The 
strong  variation  of  8  with  electric  field  ensures  that  breakdown  occurs 
near  the  pulse  center  at  the  breakdown  threshold. 

In  Fig.  6,  the  calculated  temperature  dependence  of  the  breakdown 
field  is  shown  for  sodium  chloride  at  1.06  and  10.6  pm.  The  data  re¬ 
ported  by  the  Lebedev  group  are  superimposed.  At  1.06  pm,  the  agree¬ 
ment  between  measured  and  calculated  values  is  excellent  for  variations 
in  both  magnitude  and  temperature  of  the  breakdown  field.  At  the 
lowest  temperature,  where  the  data  fall  well  below  our  calculated 
values,  the  measured  breakdown  field  may  not  be  intrinsic. 

At  10.6  pm,  the  calculated  breakdown  field  is  nearly  independent  of 
temperature,  in  accordance  with  the  data.  However,  the  calculated 
value  of  the  breakdown  field  is  smaller  than  the  reported  values  by 
a  factor  of  two.  The  smaller  value  of  the  theoretical  threshold  is 
in  sharp,  welcome  contrast  to  previous  theoretical  results,  which 
have  been  unable  to  explain  the  breakdown  at  the  low  experimental 
values  of  Eg.  More  complete  information  on  the  magnitude  and  energy 
variations  of  the  acoustic  phonon  contribution  to  y^  and  y^--particu- 
larly  their  dependence  on  Q(£j)--may  remove  that  quantitative  dis¬ 
crepancy.  The  value  of  the  breakdown  field  at  1.06  pm  is  fairly  in- 
sensitive  to  the  value  of  Q(5j)  because  Eg  is  proportional  to  the 
ratio  y^/y^,  as  discussed  in  Sec.  VI.  In  contrast,  agreement  at 
10.6  pm  is  improved  markedly  by  an  increased  value  of  Q,  because 

a 

9 

Eg  is  proportional  to  the  product  y^.  Although  it  is  feckless 
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to  adjust  the  values  of  parameters  to  force  theoretical  results  to  fit 
experimental  data,  selecting  the  value  of  a  single  parameter  can  never¬ 
theless  improve  the  agreement  between  theory  and  experiment,  as  shown 
in  Fig.  6,  for  which  an  adjusted  value  of  Q,  was  used,  rather  than  the 

a 

2 

unadjusted  value  of  Q,  =  0.35  nm  (which  is  not  accurately  known). 

a 

The  breakdown  fields  at  1.06  pm  for  several  alkali  halides  are 
summarized  in  Fig.  7.  The  present  theory  provides  an  excellent  account 
of  the  variation  of  the  breakdown  field  from  material  to  material. 

In  Fig.  8,  the  variation  of  the  calculated  breakdown  field  with 

pulse  duration  is  compared  with  the  experimental  results  of  the  Lebedev 

12  13  9 

group;  ’  Fradin,  Yablonovitch,  and  Bass;  Smith,  Bechtel,  and 

Bloembergen;^  and  Fradin,  Bloembergen,  and  Letellier.^  The  theo¬ 
retical  values  track  the  trends  in  the  data  nicely,  though  our  calcu¬ 
lated  breakdown  fields  are  somewhat  larger  than  those  measured  by  the 
Bloembergen  group. 

The  theory  adequately  accounts  for  the  data  on  magnitude,  temper¬ 
ature  dependence,  pulse-duration  dependence,  frequency  dependence,  and 
material-to-material  variation  with  no  adjusted  parameters.  There  are 
minor  qualitative  discrepancies  at  10.6  pm,  which  have  at  least  four 
possible  causes.  First,  even  though  the  description  of  yk  and  ^phYL 
involves  no  adjustable  parameters,  the  theory  is  in  fact  phenomeno¬ 
logical.  A  more  quantitative  account  of  these  basic  relaxation  rates 
may  improve  the  agreement  between  theory  and  experiment.  In  particu¬ 
lar,  the  calculated  values  are  inherently  less  accurate  at  10.6  pm  than 
at  1.06  pm  and  other  short  wavelengths--because  small -wave- vector 
phonons  determine  the  breakdown  at  1.06  pm,  whereas  phonons  with  wave 
vectors  q  s  2kg^  determine  the  breakdown  at  10.6  pm. 
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Section  VI  shows  that  low-energy  electrons  determine  the  break¬ 
down  at  1.06  pm,  whereas  high-energy  electrons  govern  it  at  10.6  ym. 

The  maximum  phonon  wave  vector  is  q  =  2k  (where  k  ~  L  is  the  elec¬ 
tron  wave  vector),  as  seen  in  the  upper  limit  of  the  integral  in  Eq. 
(4.9).  For  phonon  wave  vectors  q  near  2kg^,  the  results  become  more 
model -dependent.  In  the  calculation  of  y^  and  y^  in  Sec.  IV,  special 
treatment  is  required  to  remove  nonphysical  divergences  related  to 
Qg  and  factors  of  l/i^,  which  diverges  as  co^  approaches  zero  at  q  =  2kg^, 
A  better  treatment  of  the  large-wave-vector  phonon  scattering  should 
increase  the  theoretical  values  of  Eg  at  10.6  ym. 

Second,  the  inaccuracy  of  the  value  of  Q(£j)  discussed  above  may 

account  for  the  small  discrepancies.  Third,  low  theoretical  10.6  ym 

values  of  Eg  may  result  from  high  experimental  values  due  to  spherical 
26 

aberration,  as  discussed  in  Sec.  VI.  Fourth,  the  breakdown  at  10.6  ym 
may  be  limited  by  the  generation  of  starting  electrons,  rather  than  by 
the  avalanche  process  itself. 

Other,  less  likely  explanations  include  Raman  scattering  and  dif¬ 
fusion  of  the  electrons  out  of  the  focal  volume  (which  is  unlikely  be¬ 
cause  the  focal  volumes  at  10.6  ym  are  greater  than  those  at  1.06  ym 
and  because  simple  estimates  indicate  that  real-space  diffusion  is 
negligible).  Although  certain  low  experimental  values  of  Eg  agree 
well  with  theory,  those  that  are  lower  than  the  greatest  reliable 
value  are  currently  believed  to  be  extrinsic.  Finally,  it  is  never 
certain  that  the  measured  breakdown  field  is  the  intrinsic  value. 
However,  in  experiments  where  the  intrinsic  limit  is  not  reached, 
the  experimental  values  will  usually  be  lower  than  the  theoretical. 
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VI.  INTERPRETATIONS  AND  CALCULATIONS  USING  SIMPLE  MODELS 


Physical  explanations  of  the  results  presented  in  Sec.  V  are 
informative,  important,  and  invaluable  in  establishing  intuition.  The 
simplest  model--the  average-electron  model  discussed  in  Sec.  II— 
affords  a  useful  framework  for  discussion  of  the  results,  but  by  itself 
fails  to  account  for  experimental  observations.  Calculations  based 
on  the  transport-equation  approach  discussed  in  Secs.  Ill  and  IV  agree 
well  with  experiment,  but  numerical  solutions  to  differential  equations 
without  verification  and  interpretation  are  vacuous. 

Constant-V  and  -D  model .  This  model  treats  the  transport  equa¬ 
tion  in  Sec.  Ill  exactly  for  the  case  in  which  V  and  D  are  constants, 
independent  of  energy.  Proper  application  of  the  constant-V  and  -D 
model,  for  which  Appendix  A  gives  the  exact  solution  of  the  diffusion 
equation,  allows  the  calculation  of  Eg  from  simple  closed-form  expres¬ 
sions  that  require  no  numerical  computation.  The  resulting  accuracy  of 
these  highly  simplified  calculations  is  well  within  the  range  of 
accuracy  of  the  model  and  values  of  the  input  parameters.  Thus,  the 
model  results  afford  a  simple  method  of  obtaining  values  of  Eg,  as  well 
as  providing  scalings  and  aiding  intuition. 

The  closed-form  results  of  the  constant-V  and  -D  model  also  illus¬ 
trate  the  lowering  of  Eg  below  the  average-electron  value  in  Eq.  (2.7), 
as  well  as  other  features  of  the  general  theory.  Simple  limiting  forms 
of  the  breakdown  field  Eg  and  the  multiplication  rate  8  for  the 
constant-V  and  -D  model  derived  in  Appendix  A  are 
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EB  EdEVO  ’ 

where 

fs  [1  +  C\  g  +  (1  +  £a  g)]  1//2  ,  for  E2  «  E2Q  , 
ss  (0.219&j/g£  +  0.57)1/2  ,  for  E2  »  E2Q  , 

with 

9  5  %hYl.y51*I  <6-3) 

and 

3  fs  (e^^Y^/^aE2)  exp  (-E2Q/E2)  ,  for  E2  «  E2Q  ,  (6.4a) 

4.6/io  hy.  r\  n  99 

» - (E  '  0.57EJ0)  ,  for  l  »  E}0  .  (6.4b) 

*IEV0 

Here  Eyg  defined  in  Eq.  (2.6)  is  the  value  of  E  at  which  the  energy- 
space  velocity  is  equal  to  zero. 

Calculation  of  Eg  using  the  constant-V  and  -D  model.  In  order 
to  apply  the  constant-V  and  -D  results  to  real  problems--in  which 
and  yL  are  functions  of  energy,  rather  than  constants- 
must  be  evaluated  at  appropriate  values  of  &,  and  an  effective  value 
^Ieff  ^1  ^at  1S  srna^er  than  the  rea^  must  usually  be  chosen. 

The  magnitude  of  the  error  in  the  resulting  value  of  Eg  depends  on  how 
well  these  values  of  y^,  y^,  and  £j  (or,  equivalently,  Eyg,  yl,  and  £j) 
can  be  chosen.  Fortunately,  the  shapes  of  the  Yk(£)  and  yl(&)  curves 
are  such  that  choice  of  the  constants  EyQ,  yL,  and  are  obvious  and 


-yk(£)  and  yL(&) 


(6.1) 

(6.2a) 

(6.2b) 


the  accuracy  of  Eg  is  extremely  good,  as  illustrated  by  the  following 
examples. 

2  2 

In  the  discussion  of  V  and  D  below,  the  quantity  Eyg  -  E  is  a 
measure  of  the  height  of  the  barrier  in  energy  space  through  which  the 
electrons  must  diffuse.  Thus,  the  high  barrier  extending  from  near¬ 
zero  energy  to  L  =-  0.3<£j  for  1.06  ym  in  Fig.  9  virtually  determines 
the  value  of  Eg--because  the  value  of  E2q(&)  -  E2  as  500  (MV/cm)2  at 

=  0.15  is  a  factor  of  five  greater  than  ~100  (MV/cm)  at  &/£ j  = 
0.45,  and  the  value  of  Eg  is  sensitive  to  the  barrier  height.  With 
Evq  =  15.3  MV/cm,  £jeff  =  0.22^,  and  yl  =  1-6  *  ln'14  s"1— all 
corresponding  to  the  barrier  shown  as  the  dashed  line  in  Fig.  9-- 
Eqs.  (6.1),  (6.2a),  and  (6.3)  give  g  =  460,  Erf  =  0.33,  and 

Ed  =  5.1  MV/cm  ,  const.  V  and  D,  1.06  ym  , 

*  (6.5) 

=  5.6  MV/cm  ,  full  diff.  eq.,  1.06  ym  . 

The  approximate  value  of  5.1  MV/cm  nearly  matches  the  value  of 
5.6  MV/cm  obtained  from  the  numerical  solution  of  the  diffusion  equa¬ 
tion.  The  corresponding  results  for  10.6  ym  for  the  narrow  and  wide 

5  4 

barriers  in  Fig.  9  are,  respectively,  g  =  1.67  *  10  and  3.72  x  10  , 

Ed  =  0.25  and  0.27,  and 

Eg  =  0.72  MV/cm  and  0.70  MV/cm  ,  const.  V  and  D,  10.6  ym  , 

(6.5a) 

=  0.68  MV/cm  ,  full  diff.  eq.,  10.6  ym  , 

again  showing  excellent  agreement.  The  two  barriers  are  two  approx¬ 
imations  to  the  actual  curves.  They  are  used  to  show  the  insensi¬ 
tivity  of  the  results  to  the  choice  of  the  barrier  in  this  case. 
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2  2 

Magnitude  of  Eg.  In  the  low-field  limit  E  E^q,  which  is  valid 
for  nanosecond  pulses,  the  breakdown  field  in  Eqs.  (6.1)  and  (6.2a)  is 
less  than  E^q,  a  typical  value  being  Eg  =  0.3EVQ.  Thus,  the  simple 
model  illustrates  the  energy-space  diffusion  effect  of  reducing  the 
value  of  Eg  below  E^q  for  nanosecond  pulses.  That  is,  E^  in  Eq.  (6.1) 
is  less  than  unity. 

2  2 

In  the  high-field  limit  E  »  Eyg,  which  is  approached  for  pico¬ 
second  pulses,  the  breakdown  field  in  Eqs.  (6.1)  and  (6.2b)  is  greater 
than  EVQ.  The  high-field  limit  is  valid  for  sufficiently  short  pulse 
durations,  because  a  decrease  in  the  pulse  duration  requires  an  in¬ 
crease  in  the  electric  field  in  order  for  the  electrons  to  be  accel¬ 
erated  in  the  shorter  time.  Conversely,  at  shorter  pulse  durations 

such  as  nanoseconds,  the  electric  field  is  small,  and  the  inequality 
2  2 

E  «  Evq  is  satisfied. 

The  model  shows  that  the  average-electron  value  of  Eg  in  Eq.  (2.7) 
affords  a  very  rough  approximation  to  Eg.  Moreover,  the  corrections 
from  the  transport-equation  approach  show  Eg  to  be  somewhat  smaller 
than  Eyg  for  nanosecond  pulses  and  somewhat  greater  for  picosecond 
pulses. 

Frequency  dependence  of  Eg.  Sparks^  pointed  out  a  quite  general 
shortcoming  of  previous  theories.  The  factor  (1  +  y^/w  )  '  in  Eq. 
(2.6)  is  commonly  used  because  it  comes  from  the  conductivity  in  that 
equation.  If  xk  were  independent  of  frequency,  as  was  well  accepted 
at  that  time,  the  resulting  theoretical  frequency  dependence  of  the 
breakdown  electric  field  Eg  would  grossly  disagree  with  experimental 
results.  The  resolution  of  that  difficulty--and  the  reason  for  the 


successful  explanation  of  the  frequency  dependence  of  Eg  by  the  present 
theory--is  that  xk  in  the  simple  average-electron  model  cannot  be 
treated  as  a  constant,  independent  of  frequency.  The  value  of  xk(<l) 
is  of  course  independent  of  the  laser  frequency.  However,  as  a  result 

of  the  energy  dependence  of  xk  and  x^,  different  effective  values  of 

xk  must  be  used  at  different  laser  frequencies,  as  illustrated  by  the 
following  example. 

At  1.06  pm,  the  inequality  Yk  «  is  valid;  whereas  at  10.6  pm, 
the  inequality  w  «  Yk  is  valid.  Thus,  Eq.  (2.6)  gives 

Evo  ss  to(3/ai)phm/e2)1^2(YL/Yk)1//2  >  at  1.06  pm  ,  (6.6) 

and 

Evo  as  ^/(i*iphm/e^) (Y|_Yk)^^  »  at  10.6  pm  .  (6.7) 

Equation  (6.6)  shows  that  the  1.06  pm  breakdown  field  is  controlled 
by  the  ratio  Y[_/Yk»  which  is  greatest  at  low  electron  energies,  as 
illustrated  in  Fig.  1.  The  energy  dependence  of  Eyg,  which  is  pro¬ 
portional  to  (Y|_/Yk)  ^  at  1-06  pm,  is  plotted  in  Fig.  9.  As  men- 
tioned  above,  Eyg(£)  -  E  is  a  measure  of  the  height  of  the  barrier 
through  which  the  electrons  must  diffuse.  The  obstacle  to  accelerating 
the  electrons  to  energy  £j  occurs  at  low  energy,  where  EVg  has  the 
greatest  value  for  the  1.06  pm  laser  frequency,  as  shown  by  Fig.  9. 

In  contrast,  at  10.6  pm.  Eg  is  controlled  by  the  product  y^Yk  ,  which 
is  greatest  at  high  electron  energies  £ss£j,  as  shown  by  Figs.  1  and 
9.  Thus,  the  obstacle  occurs  at  high  energies  £  5s  for  the  10.6  pm 
laser  frequency. 


-36- 


Since  low-energy  electrons  are  important  for  inducing  breakdown 
at  1.06  ym  and  high-energy  electrons  at  10.6  ym,  the  corresponding 
values  of  vk  should  be  used  in  evaluating  EVQ  in  Eq.  (2.6).  The  small, 
low-energy  value  of  y^  from  Fig.  1  should  be  used  at  1.06  ym;  but  the 
large,  high-energy  value  of  y^  at  10.6  ym.  (The  corresponding  values 
of  Yl  must  also  be  used,  of  course--but  y^  does  not  vary  as  strongly 
with  energy  as  does  y^.)  The  diffusion  equation  automatically  accounts 
for  the  correct  at  all  frequencies;  only  when  considering  the 
average-electron  model  must  the  appropriate  be  selected  explicitly. 

For  wavelengths  shorter  than  approximately  one  micrometer, 

»  1  is  sufficiently  well  satisfied  at  all  £  for  the  shape  of  the 
Evo(£)  curves  to  be  independent  of  u>--that  is,  EVQ(£)  ~  w.  Thus,  for 
X  ~  1  ym,  the  diffusion  transport  theory  gives  Eg  ~cu.  However,  any 
of  several  effects  is  expected  to  cause  Eg  to  be  lower  than  the 
diffusion-theory  value  at  short  wavelengths.  These  effects  include 
the  breakdown  of  the  validity  of  replacing  differences  by  differentials, 
the  specific  large-quantum  processes  discussed  by  Sparks,^  and  multi¬ 
photon  absorption. 


Dependence  of  Eg  on  tp,  £j ,  fwphYL,  and  w.  In  the  low-field  limit 
2  2 

E  «  E^q,  which  is  valid  for  nanosecond  pulses,  Eqs.  (6.1)  and  (6.2a) 
show  that  Eg  is  weakly  dependent  on  the  parameters  t  and  £j  that  appear 


in  E^,  and  is  more  strongly  dependent  on  the  parameters  and  u 

that  appear  in  Eyg.  The  weak  dependence  of  the  breakdown  field  on  the 


bandgap  (or,  to  be  precise,  on  £j)  and  on  the  laser-pulse  duration  t 
for  nanosecond  pulses  was  at  first  surprising.  However,  the  previously 


used,  average-electron  criterion  for  breakdown,  Eg  =  EVq,  also  gave  Eg 
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independent  of  <Sj  and  t  .  Furthermore,  agreement  with  experimental 

2  2 

data  is  good.  In  the  high-field  limit  E  »  Eyg,  which  is  approached 
for  picosecond  pulses,  Eqs.  (6.1)  and  (6.2b)  show  Eg  to  be  rather 
strongly  dependent  on  all  parameters--tp,  &j,  A/u)phY|_»  ancl  U)- 

Focal-volume  dependence  of  Eg.  In  laser  breakdown  of  gases,  the 
diffusion  of  the  electrons  out  of  the  focal  volume  (not  to  be  confused 
with  energy- space  diffusion)  is  an  important  effect  that  determines  the 
value  of  Eg.  However,  simple  estimates  indicate  that  the  diffusion 
of  electrons  out  of  the  focal  volume  in  solids  is  negligible  in  all 
experiments  reported  to  date.  Another  factor--the  probability  of  hav¬ 
ing  or  generating  a  starting  electron  in  the  focal  volume--could  con¬ 
ceivably  cause  Eg  to  depend  on  the  focal  volume,  as  discussed  in  Sec.  I. 

Soileau  and  coworkers  suggested  that  spot-size  scaling  could 
possibly  reconcile  the  results  of  various  analysts.  They  proposed  that 
the  value  of  Eg  should  be  a  function  of  the  focal-volume  size  because 
the  probability  of  having  a  starting  electron  present  decreases  with 
volume.  Verification  of  the  spot-size  scaling  law— through,  for  ex¬ 
ample,  experiments  with  starting  electrons  supplied  by  a  separate  source, 
or  examination  of  the  probabilistic  nature  of  breakdown  in  small  focal 
volumes  and  its  deterministic  nature  in  large  focal  volumes--is  needed 
because  the  starting  electrons  are  currently  thought  to  be  easily  gen¬ 
erated  by  lattice  imperfections  such  as  impurities,  dislocations,  or 
grain  boundaries.  Breakdown  is  thought  to  be  controlled  by  the 
avalanche  process  rather  than  the  generation  of  seed  electrons.  How¬ 
ever,  these  hypotheses  have  not  yet  been  supported  either  by  experiment 
or  convincing  theory,  so  verification  of  the  suggested  mechanism  giving 
the  spot-size  scaling  would  be  useful. 
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Aaron,  Ireland,  and  Grey  Morgan  showed  that  "entirely  spurious 

'focal-length'  dependence  of  1^,  which  is  not  infrequently  found  in 

published  literature"  can  result  from  irradiance  reduction  caused  by 

spherical  aberration.  Extreme  care  in  measurement  is  required  to  avoid 

the  spurious  focal-volume  effects. 

Shape  of  the  8(E)  curves.  The  central  features  of  the  electric- 

field  dependence  of  6  illustrated  by  Figs.  4  and  5  result  from  quite 

general  considerations  that  are  independent  of  specific  models.  At 

high  fields,  the  diffusion  in  energy  space  becomes  negligible,  and  the 

electrons  stream  freely  to  the  multiplication  energy  according  to 

2 

the  relation  do/dt  (d£/dt)^  ~  oE  .  As  discussed  in  Sec.  II,  that 

2  2 

streaming  gives  8  ~  E  at  high  fields.  The  E  dependence  of  8  at  high 
fields  is  also  explicitly  illustrated  by  the  constant-V  and  -D  result 
in  Eq.  (6.5). 

At  low  fields,  8  decreases  rapidly  with  decreasing  E  because  the 
2  2 

barrier  height  Eyg(£)  -  E  through  which  the  electrons  must  diffuse 

becomes  high.  Indeed,  in  the  average-electron  model  of  Sec.  II,  8  goes 

to  zero  at  E  =  E^q.  For  the  diffusion  model,  8  goes  to  zero  rapidly 
2  2 

as  8  ~  exp  (-Ec/E  ),  as  shown  quite  generally  in  Appendix  B.  The  ex¬ 
plicit  result  in  Eq.  (6.4a)  illustrates  that  general  high-field  re¬ 
sult,  with  the  characteristic  electric  field  Ec  equal  to  Eyg. 

Temperature  dependence  of  the  8(E)  curves.  In  the  important  range 
of  electron  energies  discussed  below,  the  inequality  u>  »  Yk(£)  is 
satisfied  at  1.06  ym.  At  high  temperatures,  Yk(&)  ~T  and  /&)phYL  ~  T® 

if  the  thermal  expansion  corrections  are  ignored.  The  electron  gains 

2  2  2 

energy  from  the  electric  field  at  the  rate  2e  y^E  /m*u)  ,  according  to 
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Eq.  (2.4);  thus,  as  the  temperature  increases,  the  electron  is  accel¬ 
erated  more  efficiently  by  the  field  because  Yk  increases.  Since  the 
loss  is  nominally  temperature- independent,  B  increases  with 

temperature  for  fixed  electric  fields.  Consequently,  the  electric 
field  required  for  breakdown  decreases  as  the  temperature  is  raised, 
in  agreement  with  the  data. 

Figure  4  shows  the  variation  of  B  with  temperature  and  wavelength 
at  10.6  pm.  The  variation  of  B  with  electric  field  at  fixed  tempera¬ 
tures  is  qualitatively  similar  to  the  1.06  pm  behavior,  as  expected, 
but  its  variation  with  temperature  is  more  complex.  As  the  temperature 
increases  from  77  K  to  room  temperature,  $  increases  much  as  it  does 
at  1.06  pm.  Over  most  of  that  range,  to  >  Yk(£)  is  satisfied  for  the 
important  range  of  energies,  and  the  behavior  of  B  with  temperature 
may  be  understood  by  the  arguments  of  the  previous  paragraph.  At  room 
temperature  and  above,  the  inequality  u  <  Yk(&)  is  satisfied;  the 
electron  thus  gains  energy  from  the  electric  field  at  the  rate 
2e  E  /3m*yk,  according  to  Eq.  (2.4).  As  the  temperature  increases, 
with  Yk  nominally  proportional  to  temperature,  the  electron  is  accel¬ 
erated  less  efficiently,  and  B  decreases  with  increasing  temperature, 
in  contrast  to  its  behavior  at  1.06  pm.  Thus,  B  first  increases  with 
temperature,  reaches  a  maximum,  and  then  falls.  Since  the  criterion 
for  breakdown  outlined  below  places  the  breakdown  field  near  the  maxi¬ 
mum  at  room  temperature,  the  breakdown  field  is  fairly  independent  of 
temperature,  again  in  accordance  with  the  data. 

Energy-space  velocity  and  diffusion.  The  energy-space  velocity 
V,  which  appears  in  Eq.  (3.2)  and  is  defined  for  the  general  case  in 


Eq.  (3.3),  determines  the  rate  at  which  the  electrons  drift  to  higher 
energies  for  positive  V,  or  to  lower  energies  for  negative  V,  as  a 
result  of  the  combined  effects  of  E  acceleration  and  Y|_  loss.  The 
energy-space  diffusion  coefficient  D,  which  appears  in  Eq.  (3.2)  and 
is  defined  in  Eq.  (3.3),  determines  the  rate  at  which  electrons  dif¬ 
fuse  from  energy-space  regions  of  high  electron  concentrations  to 
neighboring  low-concentration  regions. 

The  average-electron  model  of  Sec.  II  is  equivalent  to  neglecting 

the  energy-space  diffusion.  In  extending  the  average-electron  model 

2  2 

to  include  energy-space  diffusion,  the  quantity  EVQ(£)  -  E  --when 
positive--is  a  measure  of  the  height  of  the  barrier  through  which  elec¬ 
trons  must  diffuse  to  reach  the  multiplication  energy  Ly  When 

E vo(£)  -  E  is  negative,  the  electrons  drift  to  higher  energies  even 

2  2 

in  the  absence  of  diffusion.  Neglecting  diffusion  yields  EVQ(£)  -  E  s  0 
as  the  criterion  for  breakdown,  as  indicated  in  Eq.  (2.7). 

The  value  of  EyQ(£)  -  E  at  which  breakdown  occurs  is  positive 
over  at  least  part  of  the  energy  range  0  $  ^  £j--which  indicates 
that  electron  diffusion  in  energy  space  effectively  increases  the  elec¬ 
tron  energy  to  ty  The  electric  field  thus  need  not  be  as  large  as  the 
value  in  Eq.  (2.7)  predicted  by  the  average-electron  theory.  That  dis¬ 
crepancy  is  one  reason  previous  average-electron  theories  predicted 
values  of  Eg  greater  than  measured  values. 

Physically,  the  lowering  of  Eg  by  diffusion  occurs  as  follows: 

The  process  by  which  the  electron  energy  is  increased  to  the  multipli¬ 
cation  value  involves  many  collisions  of  the  electron  with  phonons. 
Some  collisions  (yj_)  cause  the  electron  to  lose  energy,  and  some  (y^) 
to  gain  energy  from  the  electric  field. 
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The  average  effect  of  the  collisions  determines  the  energy-space 

velocity.  The  deviations  from  the  average  determine  the  energy-space 

diffusion  coefficient.  For  example,  an  electron  that  meets  with  fewer 

than  the  average  number  of  energy- loss  collisions,  together  with  an 

optimum  distribution  of  momentum-reversing  collisions—gains  energy 

faster  than  does  an  average  electron.  At  dc,  the  optimum  distribution 

of  momentum-reversing  collisions  is  none  at  all;  electrons  undergoing 

no  collisions  are  continuously  accelerated  by  the  electric  field.  For 

id  0,  the  optimum  distribution  is  a  reversal  of  the  electron  momentum 

every  cycle  of  the  electric  field  (at  the  null  in  E).  Electrons  so 

3  4  5 

accelerated  faster  than  average  have  been  called  lucky  electrons.  ’  * 
The  diffusion  equation  accounts  for  lucky  electrons  automatically,  of 


course. 
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Appendix  A 

CALCULATION  OF  THE  BREAKDOWN  FIELD  FOR  CONSTANT  D  AND  V 


Here  we  solve  the  eigenvalue  equation  obtained  from  the  diffusion 
equation  exactly  for  the  case  of  constant  diffusion  coefficient  D  and 
velocity  V.  The  solution  for  this  simple  model  affords  approximate 
values  of  the  breakdown  field  Eg  and  the  electron  multiplication  rate 
B  in  key  limiting  cases,  and  determines  their  dependence  on  the  var¬ 
iables.  The  accuracy  of  these  highly  simplified  calculations  is  con¬ 
sistent  with  that  of  the  overall  model. 

For  constant  D  and  V,  Eq.  (3.2)  becomes 

J(£,  t)  =  Vn(£,  t)  -  D  3n(£t)/3£  .  (A.l) 

Substituting  into  Eq.  (3.1),  the  diffusion  equation  in  energy  space 
becomes 

3n/3t  =  D  32n/3£2  -  V  3n/3£  .  (A. 2) 

Setting  3n/3t  =  Bn  under  the  approximation  of  Eq.  (3.6)  yields  the 
eigenvalue  equation 

D  d2n/d£2  -  V  dn/d£  -  Bn  =  0  .  (A. 3) 

Equation  (A. 3)  is  cast  into  dimensionless  form 

d2nd/ d^d  +  Ud  dnd/d£d  -  Bdnd  =  0  (A. 4) 

2 

by  dividing  t  by  xd.f  =  £j/D  and  £  by  £j  and  introducing  the  dimension¬ 
less  variables 
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nd  -  n<4/j  >  <«d  w/oj  ,  8d  ^dif  ~  ^9^d)  *  (A. 5a) 

Ed  =  E/EV0  *  Ud  =  -  WV1  -  Ed>/Ed  *  <A-5b> 

g  =  (^pbVp75!^)  •  {A-5c) 

Substituting  the  trial  solution  nd  =  exp  (ot&d)  into  Eq.  (A. 4)  and 
solving  for  a  yields 

1  / 1  ?  \i/2 

a+  =  j  Ud  ±  r  ,  where  r  =  Ud  +  8d)  .  (A. 6) 


The  boundary  conditions  for  the  case  of  constant  D  and  V  are 


n(£j)  =  0  , 

-D  dn(0)/d£  +  Vn(0)  «  -2D  dn(£j)/d£  , 


(A. 7) 


or,  in  dimensionless  form, 


nd(l)  =  0  , 

-dnd(0)/dSd  -  Udnd(0)  =  -2  dnd(l )/d^d  . 


(A. 7a) 


The  general  solution  to  Eq.  (A. 4)  that  satisfies  Eq.  (A. 7a)  is 

nd  =  exp  [a+(£d  -  1)]  -  exp  [a_(£d  -  1)]  .  (A. 8) 

Substituting  Eq.  (A. 8)  into  Eq.  (A. 7a)  and  rearranging  terms  yields 
the  transcendental  equation  for  8 


2e 


-Ud/2 


(A. 9) 


=  cosh  r  -  (Ud/2r)  sinh  r  , 
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which  reduces  to 

cosh  C(B/D)1/2£j]  =  2  ,  for  V  =  0  .  (A. 10) 

In  the  limit  1  «  exp  ^  U^,  writing  cosh  r  and  sinh  r  as 
\  [exp  (r)  ±  exp  (-r ) ] »  with  Ud/2r  =»  1  -  20d/U2  in  Eq.  (A. 9),  we 
obtain  the  limiting  cases  of 

0d  “*  U^e  ,  for  1  «  exp  ^  Ud  ,  (A.  11a) 

-4  ‘1/Ed  2 

-  Ed  e  a  ,  for  EJ  «  1  .  (A. lib) 

Substituting  Ud  =  (1  -  E^)/E^  and  0d  =  (gE^)-^  into  Eq.  (A. 11a)  and 
formally  solving  for  Ed  yields 

i  ?  2  2  r1/2 

Ed~{l  +&zg  +  *»[(1  *  E2)  /E2]} 

-  P  +  £n  g  +  Bn  (1  +  Bn  g)]"1/z  ,  g  »  1  ,  (A. 12) 

2  2 

or  E  «  Evq.  The  breakdown  field  Eg  can  then  be  obtained  using 
EB  =  EdEV0* 

The  value  of  0  in  this  small-field  limit  is  obtained  from  Eq. 

(A. lib)  by  using  the  definitions  of  0d  and  Ed  in  Eq.  (A. 5),  with 
EV0  =  <aVL/‘)1/2-  ^1S  ®lves  resu^t  of  Efl*  (6.4a) : 

0  as  (6/wphYL^2/^aE2)  exp  (-E2Q/E2)  , 


for  E2  «  E2q  . 
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In  the  other  limit  of  high  fields,  »  1  (E2  »  E2Q),  Eq.  (A. 5b) 
gives  a*  -1.  By  substituting  =  -1  into  Eq.  (A. 9)  and  expanding 
cosh  r  and  sinh  r,  the  zeroth-order  approximation  rQ  *  1.59  to  r  is 
obtained.  Rearranging  the  definition  of  r  in  Eq.  (A. 6)  gives 

6d  =  r2  -  lUd  *  (A. 13a) 

or,  as  a  zeroth-order  approximation, 

ed  ~  Bd0  =  r0  "  i  =  2-28  (A. 13b) 

using 

6  -  B^D/S2  .  2BdfephYLSE2^Ej0  .  (A.14a) 

Thus  Eq.  (A. 13b)  yields 

B0  =  4.57/fcphYL£E2/£2E20  .  (A. 14b) 

By  setting  Bq  =  51 /tp,  the  breakdown  field  can  be  calculated  from 
Eq.  (A. 14b): 


(0.219£I/g&)1/2EV0  . 


(A. 15) 


A  better  approximation  is  obtained  by  calculating  the  next-order 
term  in  B  and  Eg  using  a  double  Taylor  expansion  about  and  r  in 
Eq.  (A. 9): 


=  cosh  (rn  +  6) 


<Udn 

2(rn 


+  €) 


sinh  (rn  +  6) 


(A. 16) 


Expanding  and  equating  first-order  terms  gives  6  =  -0.336C  and 


ed  ■  <ro  +  s)2  -  T  ■  ro  -  \  *  2r06  +  7  «  ■  Bd0  +  6, 

edi  s  2ro{  +  1  e  '  ' 

Since  Ud  .  (1  -  E2)/Ejj  ■  -  1  +  1/sjj  -  Ud0  ♦  «dl,  then 

$dl  =  -  0.570/E2  =  -  0.570E2Q/E2  . 

2  2 

Thus,  in  the  limit  of  E  »  EVq, 

B-4.57^)phYL^(E2  -  0.57E20)/£2E20 
and 

Eg  -  (0.219SjM)1/2EV0  . 


. 570€  . 


-48- 


Appendix  B 

AN  ANALYTIC  EXPRESSION  FOR  THE  MULTIPLICATION  RATE  £ 
IN  THE  LIMIT  OF  A  WEAK  ELECTRIC  FIELD 


Here  we  derive  an  expression  for  the  electron  multiplication  rate  B 
appearing  in  Eq.  (3.6),  in  the  limit  that  the  laser  field  E  is  weak. 

We  begin  by  substituting  Eq.  (2.6)  into  Eq.  (3.1)  to  find 

Bn(S)  +  SJ/a5  =  0  ,  (B.l) 


which,  upon  integrating  over  energy  from  &  -  0  to  t  -  £j,  yields  the 
following  general  expression  for  £: 


B  = 


d£n(£) 


(B.2) 


We  used  the  boundary  condition  in  Eq.  (3.5)  to  obtain  this  result. 

Now  in  the  limit  of  very  small  laser  fields  E,  the  electron  multi¬ 
plication  rate  B  will  also  be  very  small--as  demonstrated  by  physical 
considerations  as  well  as  the  numerical  calculations  presented  in  Secs. 
Ill  through  V.  For  very  small  values  of  £,  the  term  Bn(£)  in  Eq.  (B.l) 
may  be  ignored  everywhere,  and--for  all  energies  save  those  near  the 
bottom  of  the  band  where  electrons  are  injected--we  have 

3J/9&  =  0  ,  (B.3) 


or 
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J(<S)  =  V(<L)n(£)  -  D(£)  9n/96  =  J(<Lj)  ,  (B.4) 


a  relation  that  holds  for  all  values  of  t  except  those  near  zero. 

Equation  (B.4)  is  an  elementary  differential  equation,  and  may  be 
integrated  at  once,  subject  to  the  boundary  condition  n (<t j )  =  0.  The 
result  is 


r  1 

r1 

V 

n(£)  =  J(£j)  exp 

f  d£" 

f 

d£* 

- exp 

-  f  d£"  ® 

mJ0  D(5")J 

l 

D(£*) 

.  j 0  . 

(B.5) 


Hereafter  we  use  Eq.  (B.5)  in  combination  with  the  low-field  forms 


m,)  =  -  fa  hyL 


(B.6a) 


and 


m  = 


_  2e2E2 


3m*  2  ,  2tv\ 

W  +  Yk(i) 


(B.6b) 


The  central  quantity  is  the  integral 


q(5)  =  f  )/D(&' )  , 

•'O 


(B.7) 


which  is  a  negative  definite  function  of  energy,  with  maximum  absolute 
value  at  6  =  £j.  In  fact,  in  very  low  fields,  we  have  |q(£j)|  »  1, 
and  also  |£jq'(£j)|  »  1.  The  second  integral  in  Eq.  (B.5)  is  con¬ 
trolled  by  the  contributions  from  energies  very  close  to  £j,  so  we 
have 
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<6 


d£* 

D(i') 


exp  [ -q (o 1 )] 


exp  [-q(ij)] 
D(^) 


exp  [ ~{L ' 


^I)q‘(wI)] 


exp  [-q(£j)]  -  exp  [ — q ) ] 

“  |v(Sj)'| 

For  n(£),  therefore,  we  have 


n(£)  = 


J(£t)  ,  i 

- —  {exp  [q(£)  -  q(£j)3  -  l} 

|v  £j)  1  1  ’ 


(B.8) 


(B.9a) 


J(£t) 

a* -  exp  [q(g)  -  q(£,)3  ,  (B.9b) 

|v(Sj)|  1 

where  the  factor  of  unity  in  £q.  (B.9a)  is  important  only  very  close 
to  £  =  £j,  a  region  that  contributes  little  to  the  denominator  in 
Eq.  (B.2). 

If  we  use  the  form  for  Eq.  (B.9b)  in  Eq.  (B.2),  then  we  have 


e  =  |v(£j)| 


exp  [q(£j)3 


d£  exp  [q(£)3 


(B. 10) 


or,  using  explicit  forms  for  the  various  quantities  given  above. 
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Table  I.  Influence  on  the  breakdown  threshold 

effects  in  the  theory. 

Eg  of  including  various 

Including  this  effect  in  the  theory  causes 

the  theoretical  Eg  to: 

Increase 

Decrease 

Spherical  aberrations 

X 

Self-focusing 

X 

Macroscopic  inclusions 

X 

Sample  surface  effects 

X 

Impurity  levels  in  the  bandgap 

X 

Control  of  Eg  by  starting  electrons 

X 

Electron  or  thermal  diffusion  out  of 

the  focal  volume 

X 

Raman  scattering 

Xa 

Xb 

Large  quantum  nature  of  fra 

X 

Multiphoton  absorption 

X 

Inter-conduction-band  transitions 

X 

NOTE:  Experimentally  measured  extrinsic  values  are  lower  than  the 
intrinsic  theoretical  values. 

aFor  increasing  focal  diameter. 

bAs  photon  frequency  increases. 
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Table  II.  Parameters  used  in  numerical  calculations. 


Parameter  Definition  Value  for  NaCl  of  293  K 

(1/31)  eV 


/<Wph  Phonon  energy 

Exciton  threshold 
ew  High-frequency 

dielectric 
constant 

Cq  Static  dielectric 

constant 

n.  Atomic  density 

a 

Halide  atomic 

d 

cross  section 

M(£)  Lattice  mass 

&P2  Electron  energy  at 

Brillouin  zone  edge 


7.8  eV 
2.34 

5.62 

-23  -3 

2.24  x  10  “  cm  0 

2 

0.35  nm 

S9.77  x  10"23  g,  l  <  \  lu 

5.93  x  10‘23  g,  £  >  \  SBZ 
4.56  eV 


m*(S) 


Electron  effective 
mass 


^  m,  te  £  2 


BZ 


m(£/SB7),  7  &n7 


'BZ' ’  2  WBZ 
m,  l  >  ^BZ 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure  6 


Figure  7 
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FIGURE  CAPTIONS 


Variation  of  the  electron-momentum  relaxation  rate  Yk  with 
energy  (solid  curve)  for  sodium  chloride  at  room  tempera¬ 
ture.  The  individual  contributions  from  Frohlich  scatter¬ 
ing  (longitudinal-optical)  and  from  scattering  by  acoustical 
phonons  are  shown  as  dashed  lines. 

Variation  of  the  electron-energy  loss  rate  Yk  with  energy 
(solid  curve)  for  sodium  chloride  at  room  temperature. 

The  individual  contributions  from  Frohlich  scattering 
(longitudinal -optical)  and  from  scattering  by  acoustical 
phonons  are  shown  as  dashed  lines. 

Variation  of  both  the  electron-momentum  relaxation  rate 
Yk  and  the  electron-energy  loss  rate  yl  with  energy  for 
sodium  chloride  at  room  temperature. 

Variation  of  the  electron  multiplication  rate  B  with  elec¬ 
tric  field  for  sodium  chloride  at  several  temperatures  at 
1.06  ym. 

Variation  of  the  electron  multiplication  rate  B  with  electric 
field  for  sodium  chloride  at  several  temperatures  at 
10.6  ym. 

1 3 

Comparison  of  the  experimental  and  theoretical  values  for 
the  temperature  dependence  of  Eg  for  sodium  chloride  at 
1 .06  ym  and  10.6  ym. 

1 3 

Comparison  of  the  experimental  and  theoretical  values  for 


the  room- temperature  breakdown  field  for  several  alkali 
halides  at  1.06  ym. 


•  13 

Figure  8.  Comparison  of  the  experimental  and  theoretical  values 
for  the  room- temperature  breakdown  field  for  sodium 
chloride  at  1.06  pm  for  various  pulse  durations. 

Figure  9.  Plot  of  E^q,  defined  in  Eq.  (2.6),  for  sodium  chloride 

at  room  temperature.  Dotted  lines  indicate  the  values  of 
the  breakdown  field  Eg  calculated  from  the  rectangular 
barriers  shown  by  dashed  lines. 
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